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We present a measurement of the direct CP- violating charge asymmetry in B mesons decaying 
to J/ipK and J/i/m where J/ip decays to /i + fi~ , using the full Run II data set of 10.4 fb _1 of 
proton-antiproton collisions collected using the DO detector at the Fermilab Tevatron Collider. A dif- 
ference in the yield of B~ and B + mesons in these decays is found by fitting to the difference between 
their reconstructed invariant mass distributions resulting in asymmetries of A J ' ^ K — [0.59 ± 0.37] %, 
which is the most precise measurement to date, and A J '^' rr = [—4.2 ± 4.5] %. Both measurements 
are consistent with standard model predictions. 
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Currently all measurements of CP violation, either in 
decay, mixing, or in the interference between the two, 
have been consistent with the presence of a single phase in 
the CKM matrix. The standard model predicts that for 
b —> sec decays, the tree and penguin contributions have 
the same weak phase, and thus no direct CP violation is 
expected in the decays of B ± mesons to J/4'K ± . Esti- 
mates of the effect of penguin loops [jj show that there 
could be a small amount of direct CP violation of up 
to O(0.3%). A measurement of a relatively large charge 
asymmetry would indicate the existence of physics be- 
yond the standard model [l|-l3j . In the transition b — > dec, 
the tree and penguin contributions have different phases, 
and there may be measurable levels of CP violation in 
the decay B ± -> J/yj^ 0,@. 

The CP-violating charge asymmetry in the decays 
S ± — >• J/i/jK and S ± — > J '/tpir are defined as 
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Previous measurements of A J ^ K [d-llO] have been aver- 
aged by the Particle Data Group with the result A J ^ K = 
[0.1 ±0.7]% [ll|]. The most precise measurement of 
A J /^ K was made by the Belle collaboration [6j , with a to- 
tal uncertainty of 0.54%. The most precise measurement 
of A J ^ 7T was made by the LHCb collaboration 12], with 
a total uncertainty of 2.9%. The LHCb measurement is 
actually a measurement of the difference, A J ^™ — A J /^ K , 
and assumes that A J ^ K is zero. The previous measure- 



ment made by the DO Collaboration [7j has a total un- 
certainty of 0.68% for A J I'^ K and 8.5% for A J I^ using 
a data sample of 2.8 fb _1 of proton-antiproton collisions. 

This Letter presents substantially improved measure- 
ments of A J ^ K and A J ^ n using the full Tevatron Run II 
data sample with an integrated luminosity of 10.4 fb~ x . 
We assume there is no production asymmetry between 
B + and B~ mesons in proton-antiproton collisions. An 
advantage of these decay modes into J/ipX is that no 
assumptions on the CP symmetry of subsequent charm 
decays need to be made. 

These updated measurements of A J ^' K and A J ^'^ 
make use of the methods for extracting asymmetries 
used in the analyses of the time-integrated flavor-specific 
scmilcptonic charge asymmetry in the decays of neutral 
B mesons [13L Il4| . We measure the raw asymmetries 
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where N 



j/4>k- 



(Nj N , 



K-* 



is the number of reconstructed 
B~ -> J/ipK- (B+ -► J/ipK+) decays, and N J/4>7T - 
(A r j/^ 7r +) is the number of reconstructed B~ — > J/ipir~ 
(B + —¥ J/t/jir + ) decays. The charge asymmetry in B ± 
decays is then given by (neglecting any terms second- 
order or higher in the asymmetry) 



Aj/tfnr _A.J/ip-rr , a 
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where Ak is the dominant correction and is the recon- 
struction asymmetry between positively, e(K + \ a nd neg- 
atively, e(K~), charged kaons in the detector (la ]: 



i-K 



e(K+)-e(K-) 
: e(K+)+e(K~)' 



(7) 



The correction Ak is calculated using the measured kaon 
reconstruction asymmetry as described below [14|. As 
discussed later, data collected using regular reversals of 
magnet polarities results in no significant residual track 
reconstruction asymmetries, and hence, no correction for 
tracking asymmetries or pion reconstruction asymmetries 
need to be applied, hence A^ = 0. 

The DO detector has a central tracking system, con- 
sisting of a silicon microstrip tracker and the central 
fiber tracker, both located within a 2 T superconduct- 
ing solenoidal magnet ill [l6| ■ A muon system, cover- 
ing | rj | < 2 [l7|, consists of a layer of tracking detec- 
tors and scintillation trigger counters in front of 1.8 T 
toroidal magnets, followed by two similar layers after the 
toroids [18] . 

The polarities of the toroidal and solenoidal magnetic 
fields are reversed on average every two weeks so that the 
four solenoid-toroid polarity combinations are exposed to 
approximately the same integrated luminosity. This al- 
lows for a cancelation of first-order effects related to in- 
strumental asymmetries. To ensure optimal cancelation 
of the uncertainties, the events are weighted according to 
the number of J/iph^ decays for each data sample cor- 
responding to a different configuration of the magnets' 
polarities (polarity- weighting). The weighting is based 
on the number of events that pass the selection criteria 
and the likelihood selection (described below) and that 
are in the J '/iph invariant mass range used to fit the 
data. 

The data are collected with a suite of single and 
dimuon triggers. The selection and reconstruction of 
J/iph^ 1 events where h^ is any stable charged hadron and 
J/ijj — > /i + /i~ requires three tracks with at least two hits 
in both the silicon microstrip tracker and the central fiber 
tracker. The muon selection requires a transverse mo- 
mentum pt > 1.5 GeV/c with respect to the beam axis. 
One of the reconstructed muons is required to have hits 
in at least two layers of the muon system with segments 
reconstructed both inside and outside the toroid. The 
second muon is required to have hits in at least the first 
layer of the muon system. The muon track segment has 
to be matched to a particle found in the central tracking 
system. The dimuon system must have a reconstructed 
invariant mass between 2.80 and 3.35 GeV/c 2 consistent 
with the J/tp mass, 3.097 GeV/c 2 [llj]. 

An additional charged particle with px > 0.7 GeV/c is 
selected. Since the DO detector is unable to distinguish 
between K^- and 7r ± , and since the Jt\>K ± process is 
dominant, this particle is assigned the charged kaon mass 
and is required to be consistent with coming from the 



same three-dimensional vertex as the two muons, with 
the x 2 °f the vertex fit being less than 16 for 3 degrees 
of freedom. The displacement of this vertex from the 
primary proton-antiproton interaction point is required 
to exceed 3 standard deviations for the resolution of the 
vertex position in the plane perpendicular to the beam 
direction. 

The B ± selection is further improved using a likeli- 
hood ratio method taken directly from Refs. [19H22J that 
combines a number of variables to discriminate between 
signal and background: the smaller of the transverse mo- 
menta of the two muons; the x 2 of the B decay vertex; 
the B ± decay length divided by its uncertainty; the sig- 
nificance, Sb, of the reconstructed B meson impact pa- 
rameter; the transverse momentum of the h ; and the 
significance, Sk, of the h impact parameter. 

For any particle i, the significance Si is defined as 
Si = y/[e T /cr(€T)} 2 + [e L /a(e L )} 2 , where e T (ez,) is the 
projection of the impact parameter on the plane perpen- 
dicular to (along) the beam direction, and <j{zt) W{^l)\ 
is its uncertainty The trajectory of each B ± is formed 
assuming that it passes through the reconstructed B ± 
vertex and is directed along the reconstructed B mo- 
mentum. 

The final requirement on the likelihood ratio variable is 
chosen to minimize the statistical uncertainty on A-aw ■ 
The measurement of i r „ makes use of a different selec- 
tion on the likelihood ratio that minimizes the statistical 
uncertainty of ^4 r aw ■ The asymmetry results extracted 
with both of these likelihood selections are consistent. 
No event has more than one possible track and J/tp mass 
combination that passes all of the selection criteria. 

From each set of three particles fulfilling these require- 
ments, a J/tph candidate is constructed. The momenta 
of the muons are corrected by constraining the J/tp mass 
to the world average [llj . 

The number of signal candidates are extracted from the 
J/iph ± mass distribution using an unbinned maximum 
likelihood fit over a mass range of 4.98 < M(J/tjjh ± ) < 
5.76 GeV/c 2 as shown in Fig. [TJ The dominant peak 
consists of the overlap of the B^ — > J/ipK^ 1 and the 
S ± — > J/tpn^ (where the 7r ± is mis-identified as a K^) 
components. The mis-identified B ± — > J/-07r ± decay 
mode appears as a small peak shifted to a slightly higher 
mass than the B ± . The B ± —> J/ipK^ signal peak is 
modeled by two Gaussian functions constrained to have 
the same mean but, with different widths and normaliza- 
tions to model the detector's mass resolution, GK{rn). 
Taking account the DO momentum scale, the mean is 
found to be consistent with the PDG average of the B ± 
meson mass. To obtain a good fit to the data, the widths 
have a linear dependence on the kaon energy. We as- 
sume that the mass distribution of the B ± — > J/ipn^ 1 is 
identical to that of B ± — >• J/ipK^, if the correct hadron 
mass is assigned. To model the J/-07T mass distribu- 



tion, G 7r (?7i), the J/ipw signal peak is transformed by 
assigning the pion track the charged kaon mass. Partially 
reconstructed decays such as B x — ► J/tph ± X where h^ 
is any stable charged hadron and X is additional charged 
or neutral particles (e.g., the decay B ± — > J/ipK* ± ) can 
be empirically modeled with a threshold function that ex- 
tends to the B ± mass and is based on Monte Carlo simu- 
lations [20]: T(m) — arctan [pi(mc 2 — P2)] +P3, where Pi 
are fit parameters. In the default fit only the normaliza- 
tion of T(m) is allowed to vary and the other parameters 
are fixed to the values obtained from simulation. The 
combinatorial background is described by an exponential 
function, E(m), with a slope that depends on the kaon 
energy. The fractions of the J /ipK, J/ipff, and partially 
reconstructed decays depend on the h momentum. Em- 
pirical studies of the data show that this dependence can 
be modeled by the same polynomial function with dif- 
ferent scaling factors for the J/ipK, J/ipir, and partially 
reconstructed fractions. The coefficients of the polyno- 
mial and the scaling factors are determined from the fit. 
The likelihood function is defined to simultaneously fit 
the raw asymmetries, A r aw , the asymmetry of the 
partially reconstructed decays, At, and the asymmetry 
in the combinatorial background, Ae' 

C = (l- q h A J T [S K ) G K (m) + (l - q h A J x lt) G„(m) 

+ (l-q h A T )T(m) + (l-q h A E )E(m), (8) 

where q^ is the charge of the hadron. 

The raw asymmetries are extracted by fitting the re- 
sulting data sample using the unbinned maximum likeli- 
hood fit described above. The resulting J/iph^ polarity- 
weighted invariant mass distribution is shown in Fig. [T] 
The B ± -> J/tpK ± signal contains 105562 ± 370 (stat) 
events, and the B ± — > J/ipn^ 1 signal contains 3110 ± 
174 (stat) events. 

The quality of the fit is estimated by projecting the 
resulting unbinned likelihood fit onto the J/ipK invari- 
ant mass distribution (65 bins in total). A \ 2 is then 
calculated with a value of 76.2 for 47 degrees of freedom 
(the number of bins less the number of fit parameters 
excluding the asymmetry parameters). 

The invariant mass distribution of the differences, 
N(J/i(;h~)-N{J/iph + ), is shown in Fig.Hwith a result- 
ing x 2 of 58.5 for 61 degrees of freedom. The resulting 
raw asymmetries are extracted from the data are: 

A lLt K = [-0-46 ± 0.36 (stat)] %, (9) 

^r J /w * = [-4-2 ± 4.4 (stat)] %. (10) 

The background asymmetries are also determined: At — 
[-1.3 ± 1.0 (stat)] % and A E = [-1.1 ± 0.6 (stat)] %. 

The systematic uncertainties in the fitting method are 
evaluated by varying the fitting procedure. For each of 
the following variations the systematic uncertainty is as- 
signed to be half the maximum variation in the central 
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FIG. 1. The polarity-weighted J/iph invariant mass dis- 
tribution, where the h is assigned the charged kaon mass, 
after selecting on the likelihood-ratio function optimized for 
-Araw • The bottom panel shows the fit residuals (the error 
bars represent the statistical uncertainty). Fit described in 
the text. 



value. The mass range of the fit is modified so that the 
lower edge is varied from 4.95 to 5.01 GeV/c 2 , and the up- 
per edge from 5.73 to 5.79 GeV/c 2 , in 10 MeV/c 2 steps. 
This results in an uncertainty in A T l$ of 0.022% and 
in Arlt" of 0.55% (labeled "Mass range" in Table [J). 
The following modifications are made to the functions 
used to model the data. The mean of the Gaussian func- 
tions is allowed to depend linearly on the energy of the 
kaon. The pr(-?0-dependence of the width of the Gaus- 
sian function is modeled with a quadratic and a cubic 
polynomial. The parameters of the threshold function 
are allowed to float. The ratio of branching fractions for 
the decays B ± — > J/ipK^ and B ± — > J/ipn^ which are 
not constrained in the default fit are fixed to the current 
ratio from the Particle Data Group, 0.0482 [ll|, and the 
latest measurement by the LHCb experiment, 0.0381 [12J . 
This results in an uncertainty in A T l$ of 0.011% and 
in Arl^ of 0.69% (labelled "Fit function"). The effect 
of the event weighting is studied by varying the number 
of events for each magnet configuration by the statisti- 
cal uncertainty (\A/V). This results in uncertainties of 
less than 0.0005% in A J ^ K and 0.014% in A J /^, which 
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FIG. 2. The fit to the difference distribution for the data 
optimized for A T ^ (the fit is described in the text). 



are small compared to the other uncertainties and is not 
included in the summary table. 

The resulting systematic uncertainties are added in 
quadrature to obtain: 



A-J/^ 



K 



-0.46 ± 0.36 (stat) ± 0.025 (syst)] %, (11) 



Ail^ = [-4.2 ± 4.4 (stat) ± 0.88 (syst)] %. (12) 

As a cross-check the following variations of the vari- 
ous asymmetry models are also examined. The asymme- 
tries representing the threshold function and the combi- 
natoric background are set to the same value, At = Ae- 
The asymmetry of the combinatoric background is set to 
zero, Ae = 0. The asymmetry of the threshold function 
is set to zero, At — 0. The asymmetries representing 
the threshold function and the combinatoric background 
are both set to zero, Ae — At = 0. When extracting 
A-aw , the asymmetry A x l$ n is set equal to zero. When 
extracting Aaw, the asymmetry A r lw is set equal to 
zero. This results in variations in A T lw of 0.038% and 
in A r /w of 1.59%. Given the statistical and systematic 
uncertainties, the observed variations are consistent with 
no significant biases. 

To test the sensitivity of the fitting procedure, the 
charge of the charged hadron in the data is random- 
ized to produce samples with no asymmetry, and 1000 
trials are performed, each with the same statistics as the 
measurement. The central value of the asymmetry distri- 
bution, (+0.008 ± 0.011)%, is consistent with zero with 
a width of 0.37%, consistent with the statistical uncer- 
tainty found in data. These studies are repeated with 
introduced asymmetries of —1.0, —0.5 and 1.0%, each of 
which returns the expected asymmetries and statistical 
uncertainties with no significant bias. 

The residual detector tracking asymmetry has been 



j^*± _^. Kgir^ 1 decays. No significant residual track re- 
construction asymmetries are found and no correction for 
tracking asymmetries need to be applied. The tracking 
asymmetry of charged pions has been studied using MC 
simulations of the detector. The asymmetry is found to 
be less than 0.05%, which is assigned as a systematic 
uncertainty (labeled AAtracking)- 

The correction Ak (Eq. [7]), is calculated using the 
measured kaon reconstruction asymmetry presented in 
Ref. 14|. Negative kaons can interact with matter to 



produce hyperons, while there is no equivalent interac- 
tion for positive kaons. As a result, the mean path length 
for positive kaons is larger, the reconstruction efficiency 
is higher, and the kaon asymmetry, Ak, is positive. 

The kaon asymmetry is measured using a dedicated 
sample of K*°(K*°) -» K + -k~{K~-k + ) decays, based on 
the technique described in Ref. [23]. The K + tt~ and 
K~ir + signal yields are extracted by fitting the charge- 
specific Af(K ± 7r =F ) distributions, and the asymmetry is 
determined by dividing the difference by the sum. The 
track selection criteria are the same as those required for 
the J/iph^ signal. 

As expected, an overall positive kaon asymmetry of 
approximately 1% is observed. A strong dependence on 
kaon momentum and the absolute value of the pseudo- 
rapidity is found, and hence the final kaon asymmetry 
correction to be applied in Eq. [5] is determined by the 
polarity-weighted average of Ak[p(K),\t](K)\] over the 
p(K) and |t/(-?0| distributions in the signal events. A rel- 
ative systematic uncertainty of 5% is assigned to each bin 
to account for possible variations in the yield when dif- 
ferent models are used to fit the signal and backgrounds 
in the K*° mass distribution. Following studies over a 
range of fit variations, a relative systematic uncertainty 
of 3% on the J/ipK ± yields is applied. The resulting 
kaon correction is found to be (the uncertainty is labeled 
AA K in Table U): 



.4 



A' 



[1.046 ±0.043 (syst)] %. 



(13) 



studied in Ref. [13] |lj, |23| using K° s -» tt+tt" and 



The value of Ak is consistent with that presented in 
Ref. 7 taking into account the changes in the data se- 
lection and the resulting changes in the p(K) and |??(-?0| 
distributions. 

The final uncertainties are summarized in Table U 
where their combination assumes that they are uncor- 
rected. We obtain final asymmetries of 

A J/i,K = jq 5Q ± q 3 g ( gtat ) ± q Q7 ( gyst )] %) ( M ) 

A J/i>n = j_ 4 2 ± 4 4 (stat ) ± 9 ( sygt )] %> (15 ) 

This is the most precise measurement of A J ^ K to date 
and is a reduction in uncertainty by approximately a fac- 
tor of two from the previous DO result [7| . 

Several consistency checks are performed by dividing 
the data into smaller samples using additional selections 



TABLE I. The statistical and systematic uncertainties for ex- 
tracting the asymmetries A J ^ K and A J ^ 1T . 



Type of uncertainty 


A .J/i,K ( % ) 


A J/i,-n ( % ) 


Statistical 


0.36 


4.4 


Mass range 
Fit function 

^ji-t racking 

AA K 


0.022 
0.011 
0.05 
0.043 


0.55 
0.69 
0.05 

n/a 


Total systematic uncertainty 


0.07 


0.9 


Total uncertainty 


0.37 


4.5 



based on data-taking periods, magnet polarities, trans- 
verse momentum, and rapidity of the charged track rep- 
resenting the kaon. The resulting variations of A J I^ K 
and A J ^ n are statistically consistent with the results of 
Eqs.[Hand[l5l 



In summary, we have presented the most precise mea- 
surement to date of the charge asymmetry A J ^ K = 
[0.59 ± 0.36 (stat) ± 0.07 (syst)] % using 10.4 ftr 1 of 
data. In addition we have improved our measurement 
of A J I^ = [-4.2 ± 4.4 (stat) ± 0.9 (syst)] %. Both mea- 
surements are in agreement with standard model predic- 
tions. 
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